Technical  report:  Regulation  of  DNA  metabolism  by  DNA-binding  proteins  probed  by 
single  molecule  spectroscopy 


A.  Physiological  background  of  DNA  denaturation 

The  Watson-Crick  double-helix  is  the  thermodynamically  stable  configuration  of  a 
DNA  molecule  under  physiological  conditions  (normal  salt  and  room/body  temperature). 
This  stability  is  effected  (a)  by  Watson-Crick  H-bonding,  that  is  essential  for  the 
specificity  of  base  pairing,  i.e.,  for  the  key-lock  principle  according  to  which  the 
nucleotide  Adenine  exclusively  binds  to  Thymine,  and  Guanine  only  to  Cytosine.  Base¬ 
pairing  therefore  guarantees  the  high  level  of  fidelity  during  replication  and  transcription, 
(b)  The  second  contribution  to  DNA  helix  stability  comes  from  base-stacking  between 
neighboring  base  pairs:  through  hydrophobic  interactions  between  the  planar  aromatic 
bases,  that  overlap  geometrically  and  electronically,  the  base-pair  stacking  stabilizes  the 
helical  structure  against  the  repulsive  electrostatic  force  between  the  negatively  charged 
phosphate  groups  located  at  the  outside  of  the  DNA  double-strand.  While  hydrogen 
bonds  contribute  only  little  to  the  helix  stability,  the  major  support  comes  from  base¬ 
stacking. 

The  capability  of  the  two  complementary  strands  of  DNA  to  move  apart  and  to  rejoin, 
without  damaging  the  chemical  structure  of  the  two  single-strands,  is  crucial  to  many 
physiological  processes  such  as  replication  via  the  proteins  DNA  helicase  and 
polymerase,  and  transcription  through  RNA  polymerase.  During  these  processes,  the 
proteins  unzip  a  certain  region  of  the  double-strand,  to  obtain  access  to  the  genetic 
information  stored  in  the  bases  in  the  core  of  the  double-helix.  This  unzipping 
corresponds  to  breaking  the  hydrogen  bonds  between  the  base-pairs.  Classically,  the 
so-called  melting  and  reannealing  behavior  of  DNA  has  been  studied  in  solution  in  vitro 
by  increasing  the  temperature,  or  by  titration  with  acid  or  alkali  (Figure  1).  During 
thermal  melting,  the  stability  of  the  DNA  duplex  is  related  to  the  content  of  triple 
hydrogen-bonded  G-C  base-pairs:  the  larger  the  fraction  of  G-C  pairs,  the  higher  the 
required  melting  temperature  or  pH  value.  Thus,  under  thermal  melting,  dsDNA  starts  to 
unwind  in  regions  rich  in  A-T  base-pairs,  and  then  proceeds  to  regions  of  progressively 
higher  GC  content.  Conversely,  molten,  complementary  chains  of  single-stranded  DNA 
(ssDNA)  begin  to  reassociate  and  eventually  reform  the  original  double-helix  under 
incubation  at  roughly  25°  below  the  melting  temperature  . 

Complementary  to  thermal  or  pH  induced  denaturation,  ssDNA  can  be  driven  toward 
denaturation  mechanically,  by  applying  a  tensional  stress  along  the  DNA  in  an  optical 
tweezers  trap.  The  force  per  extension  increases  in  worm-like  chain  fashion,  until  a 
plateau  at  approximately  65  pN  is  reached  (Figure  2).  This  plateau  is  sometimes 
interpreted  as  new  DNA  configuration,  the  S  form.  By  a  series  of  experiments,  it  appears 
more  likely  that  the  plateau  corresponds  to  the  mechanical  denaturation  transition.  To 
first  order,  the  effect  of  the  longitudinal  pulling  with  a  force  F  translates  into  an  external 
torque  r  =  ^(F) ,  whose  effect  is  to  decrease  the  free  energy  of  melting  AG : 

^G{F)  =  ^G{F  =  0)  -  t{F)0^  (1) 

where  0^^  =  ItZ /10.35  i.S  tho  annlp  np>r  hn  nf  thp  Hni  ihip  holiv  u/lthm  it  annllpH  tnrces. 
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Figure  1.  Fraction  6^  of  double-helical  domains  within  the  DNA  as  a  function  of 
temperature  T .  Schematic  representation  of  (r) ,  showing  the  increased  formation  of 
bubbles  and  unzipping  from  the  ends  until  full  denaturation  has  been  reached. 
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model  for  a  DNA  molecule  at  the  overstretching  transition 

Figure  2.  Top:  Overstretching  transition  of  double-stranded  DNA.  The  force  extension 
relation  exhibits  a  rapid  worm-like  chain  increase  until  at  around  65  pN  a  plateau  is 
reached.  The  transition  on  the  plateau  corresponds  to  progressive  mechanical 
denaturation  from  dsDNA  (left  end)  to  ssDNA  (right  end).  Bottom:  Schematic  depiction  of 
DNA  conformation  on  the  overstretching  plateau  with  local  denaturation  zones  (bubbles). 
See  text  for  details.  Figure  courtesy  Mark  C.  Williams. 


While  the  double  helix  is  the  thermodynamically  stable  configuration  of  the  DNA 
molecule  below  the  melting  temperature  (at  nondenaturing  pH),  even  at  physiological 

conditions  there  exist  local  denaturation  zones,  so-called  DNA-bubbles,  predominantly  in 
A-T-rich  regions  of  the  genome.  Driven  by  ambient  thermal  fluctuations,  a  DNA  bubble  is 
a  dynamical  entity  whose  size  varies  by  thermally  activated  zipping  and  unzipping  of 
successive  base  pairs  at  the  two  forks  where  the  ssDNA  bubble  is  bordered  by  the 
dsDNA-helix.  This  incessant  zipping  and  unzipping  leads  to  a  random  walk  in  the 
bubble-size  coordinate,  and  to  a  finite  lifetime  of  DNA-bubbles  under  non-melting 
conditions,  as  eventually  the  bubble  closes  due  to  the  energetic  preference  for  the 
closed  state.  This  so-called  DNA  breathing  typically  opens  up  a  few  bps.  It  has  been 
demonstrated  recently  that  by  fluorescence  correlation  methods  the  fluctuations  of  DNA 
bubbles  can  be  explored  on  the  single  molecule  level,  revealing  a  multistate  kinetics  that 
corresponds  to  the  picture  of  successive  zipping  and  unzipping  of  single  base  pairs.  At 
room  temperature,  the  characteristic  closing  time  of  an  unbounded  base  pair  was  found 
to  be  in  the  range  10  to  100  psec  corresponding  to  an  overall  bubble  life  time  in  the 
range  of  a  few  msec.  The  multistate  nature  of  the  DNA-breathing  was  recently  confirmed 
by  a  UV-light  absorption  study. 

The  presence  of  fluctuating  DNA-bubbles  is  essential  to  the  understanding  of  the 
binding  of  single-stranded  DNA  binding  proteins  (SSBs)  that  selectively  bind  to  ssDNA 
and  that  play  important  roles  in  replication,  recombination  and  repair  of  DNA.  One  of  the 
key  tasks  of  SSBs  is  to  prevent  the  formation  of  secondary  structure  in  ssDNA.  From  the 
thermodynamical  point  of  view  one  would  therefore  expect  SSBs  to  be  of  an  effectively 
helix-destabilizing  nature,  and  thus  to  lower  .  However,  it  was  found  that  neither  the 

gp32  protein  from  the  T4  phage  nor  E.coli  SSBs  do.  An  explanation  to  this  apparent 
paradox  was  suggested  to  consist  in  a  kinetic  block,  i.e.,  a  kinetic  regulation  such  that 
the  rate  constant  for  the  binding  of  SSBs  is  smaller  than  the  one  for  bubble  closing.  This 
hypothesis  could  recently  be  verified  in  extensive  single  molecule  setups  using 
mechanical  overstretching  of  dsDNA  by  optical  tweezers  in  the  presence  of  T4  gene  32 
protein  [1 ,2].  Our  study  consisted  in  developing  a  biophysical  model  of  this  system,  as 
detailed  below;  see  references  [3-7],  and  references  cited  therein. 


B.  Probabilistic  modeling  by  the  master  equation  (ME) 

DNA  breathing  in  presence  of  SSBS  is  described  by  the  probability  distribution 
P[m,n,t)  to  find  at  time  t  a  bubble  with  m  open  bps  and  with  n  SSB  units  bound  to 

the  two  arches  of  the  bubble.  The  time  evolution  of  P[m,n,t)  is  given  by  a  master 
equation 

—  P(^m,n,t)  =  W  P[m,n,t)  (2) 

where  the  transfer  matrix  W  incorporates  rates  t~  for  bubble  increase  /  decrease  as 
well  as  rates  for  SSB  binding  /  unbinding: 


(3) 


The  bubble  part,  according  to  the  Poland-Scheraga  model  for  DNA-melting,  is  given  by 

^bubble {^)  =  ^0  u'”  {m  +  !)■'  ,  m>\  ,  (6) 

with  the  bubbie  initiation  factor  (j^ ,  the  weight  u  for  breaking  a  base-pair,  and  the  loop 
closure  factor  (w  + 1)  for  creating  a  polymer  loop  of  size  m  ,  with  the  offset  by  1  due 
to  persistence  length  corrections.  We  choose  typical  values,  cTp  =10'^  and  c  =  1.76 . 
Equation  (5)  is  completed  by  the  boundary  condition  (O)  =  1  •  The  contribution  from 
the  SSBs  has  the  form 

^  k"  Q.{m,n)  (7) 

with  the  binding  strength  k  =  ,  involving  the  SSB  concentration  Cq  and 

equilibrium  binding  constant  ,  where  is  the  typical  SSB  volume  and 

£'ssb  its  binding  energy.  The  weight  Q(w,n)  is  a  combinatorial  factor,  counting  aii 

possible  ways  of  putting  n  SSBs  onto  the  two  arches  of  the  bubble,  taking  into  account 
gaps  between  SSBs  due  to  the  fact  that  their  size  X  is  larger  than  a  base.  Solving  the 
above  equations  results  in  the  rates  for  the  bubble  size  f : 


t  {m,n)  =  ku 


m>l 


(8) 


t  {m,n) 


r^i+iY 

^  Q(m-l,n) 
Q.{ni,n) 


with  the  bubble  initiation  rate  t  (0,0)  -  2 ‘"ka^u.  The  SSB  number  transfer  rates  r 
result  as 


(w,«) 

r“ 


=  /C  (rt  +  l) 


Q(w,n  +  l) 
Q(w,«) 


-  nyk 


(9) 


with  the  dimensionless  ratio  y-q/k  of  the  SSB  unbinding  rate  q  and  the  base  pair 
zipping  rate  k . 

The  master  equation  (3)  can  be  solved  by  means  of  an  eigenmode  expansion  of  the 
form 


P{m,n,t)  =  Y.^pQp{m,n)  exp(r/rj 

p 


(10) 


in  which  the  coefficients  of  a  given  eigenmode  p  are  determined  via  the  initial 
conditions.  The  corresponding  eigenvalue  equation  for  the  bubble  size-SSB  number 
eigenfuction  determines  the  mode  relaxation  times  and  can  be  solved 
numerically. 


C.  Effective  binding  free  energy  in  the  limit  of  fast  SSB  binding 

Fast  SSB-binding  corresponds  to  the  limit  in  which  the  dimensionless  ratio  y  =  qlk  of 
the  SSB  unbinding  rate  q  and  the  base  pair  zipping  rate  k  is  large,  i.e.,  yU  I .  This 
limit  allows  one  to  average  out  the  SSB-dynamics  and  to  calculate  an  effective  free 
energy,  in  which  the  bubble  dynamics  with  the  slow  variable  m  runs  off.  The  result  for 
two  different  binding  strengths  k  is  shown  in  Figure  3,  along  with  the  free  energies 
corresponding  to  keeping  n  fixed.  It  is  distinct  that  while  for  lower  k  the  presence  of 
SSBs  diminishes  the  slope  of  the  effective  free  energy,  for  larger  k  the  slope  becomes 
negative.  This  implies  that  in  the  first  case  the  bubble  opening  is  more  likely,  but  still 
globally  unfavorable.  In  the  latter  case,  the  presence  of  SSBs  indeed  leads  to  full 
denaturation.  One  observes  distinct  finite  size  effects  due  to  A  >  1 :  only  when  the 
bubble  reaches  a  minimal  size  m>X,  SSB-binding  may  occur,  a  second  SSB  is 
allowed  to  bind  to  the  same  arch  only  once  m>2X ,  etc.  This  effect  also  produces  the 
nucleation  barrier  for  full  denaturation  in  the  lower  plot  of  Figure  3. 
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